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Abstract 


In this paper, we focus on the reduced basis methodology in the context of non-linear non-afhnely parametrized 
partial differential equations in which afhne decomposition necessary for the reduced basis methodology are not 
obtained ii. To deal with this issue, it is now standard to apply the EIM methodology [1,[^ before deploying 
the Reduced Basis (RB) methodology. However the computational cost is generally huge as it requires maw 
finite element solves, hence making it inefhcient, to build the EIM approximation of the non-linear terms [9|, [ij. 
We propose a simultaneous EIM Reduced basis algorithm, named SER, that provides a huge computational 
gain and requires as little as iV -|-1 finite element solves where N is the dimension of the RB approximation. The 
paper is organized as follows: we first review the EIM and RB methodologies applied to non-linear problems 
and identify the main issu^ then we present SER and some variants and finally illustrates its performances in 
a benchmark proposed in (9|. To cite this article: C. Daversin, C. Prud’homme, C. R. Acad. Sci. Paris, Ser. 
I 340 (2015). 

Resume 

Une methode EIM Bases Reduites simultanee pour les equations aux derivees partielle non- 
lineaires et non-affines. 

Dans ce papier, nous nous interessons a la methodologie bases reduites (RB) dans le contexte d’equations aux 
derivees partielles parametrisees non-lineaires et non-affines et pour lesquelles la decomposition affine necessaire 
a la methodologie RB ne pent etre obtenue id . Pour traiter ce probleme, il est a present standard d’appliquer la 
methodologie EIM dd avant de deployer la methodologie RB. Cependant le cout de calcul de cette approche est 
en general considerable car il requiert de nombreuses evaluations element fini, la rendant tres pen competitive, 
pour construire I’approximation EIM des termes non-lineaires ifl. Nous proposons I’algorithme SER qui 
construit simultanement I’approximation EIM et RB, fournit ainsi un gain de calcul considerable et requiert an 
minimum iV -|- 1 resolutions element fini ou N est la dimension de I’approximation RB. Le papier est organise 
comme suit : tout d’abord nous passons en revue les methodes EIM et RB appliquees aux problemes non- 
lineaires et identifions la difhculte principale, puis nous presentons SER et quelques variantes et finalement nous 
illustrons ses performances sur un benchmark propose par d. 

Pour citer cet article : C. Daversin, C. Prud’homme, C. R. Acad. Sci. Paris, Ser. I 340 (2015). 


1. Introduction 

Reduced order modeling is more and more used in eng ineering problems due to efficient evaluation of 
quantities of interest. The Reduced Basis Method (see EiSaiB) has been especially designed for 
real-time and many-query contexts, and cover a large range of problems among which non-afhnely parametrized 
Partial Differential Equations (PDE). A core enabler of this method is the so-called offline/online decomposition 
of the problem. This allows for computing costly parameter independent terms that depend solely on the Hnite 
element dimension. However such decomposition is not necessarily or readily available in particular for non- 
afhne/non-linear problems. The Empirical Interpolation Method (EIM, see [8|,[9j) has been developed to recover 
this core ingredient and is used prior to the reduced basis methodology on industrial based applications, see e.g. 
[I|. However the EIM building step can be costly when the terms are non-linear and requires many non-linear 
hnite element solves. It is a deterring trait for this methodology which forbids its application to non-linear 
applications. 

In this paper, we propose a Simultaneous EIM-RB (named SER) construction that requires only but a few 
hnite element solves and which builds together the afhne decomposition as well as the RB ingredient. To start we 
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first give an overview of both Reduced Basis and Empirical Interpolation methods in non-affinely parametrized 
PDE context to recall the necessary notions. Based on these considerations, a second part makes an assessment 
of the proposed simultaneous approach discussing the changes to be made in the EIM offline step. The last part 
illustrates our method with preliminar results obtained on a benchmark introduced in Q . 


2. Preliminaries 

Let u(/i) be a solution of a non-linear parametrized PDE, where /x is the p-vector of inputs G T) C MP 
the parameter space). Considering the Hilbert space X = 7Jq(D), the variational formulation of the PDE 
consists in finding u{fj,) G X as a root of a functional r such that 

fi) = 0 \fv G X. (1) 

We denote Af the dimension of the finite element approximation space X-^ C X, and the associated 

finite element approximation of u(/x) solution of ([Til . 

In the following, we use a Newton algorithm to deal with the non-linearity. Denoting ^u(/x) the solution 
at Newton’s iteration, and j the Jacobian associated with the functional r, the problem consists in finding 
— ^u{fji) G X-^ and thus G X-^ such that 

j{S’^+^u{n), V- /x; ^u{n)) = v; fi) \/v G X^ (2) 


at each iteration k. 

If j and r depend affinely in /x, there exists positive integers Qj and Qr, /x-dependent functions 9j ^, 0^'" and 
functions jn ^ (1 ^ ^ Qr, 1 ^ qj ^ Qj) allowing to write such affine decompositions 

Qj Qr 

j (m, R; ^ (u, v) r(^u, v; fi) = ^ 9'^/ {^u; (v) Vu, v G X^, M^gV (3) 

qj = l gr = l 

Eor non-afhnely parametrized problems, j and r depend on a set x; /x)}^ of non-afhne functions 

j{u, V- /x; '=m(/x); {wi{u{(i), x; /x)}i)j'=+^u(/x) = v\ /x; {wi(zi(/x), x; /x)}i) (4) 

and the affine decomposition (|3]) doesn’t exist owing to Wi(u(fj,),x', /x) definition. In order to recover an affine 
decomposition for (|3]), we build an affine approximation Wi^Mi(u{iJ,), x; /x) of x; /x) 

Mi 

Wi^MiiuifJ.),x;fi) = (i)qi^mix) Vi (5) 

m—1 


2.1. Empirical Interpolation Method 


M 

In order to build an affine approximation wm{u{pl), x; /x) = ^ [3^fj,)qm{x) of a non-affine parameter 

_ m—l 

dependent function w, we first introduce a sample Sm = {Ai> ■ • ■ > Am} Q , and the associated fonction space 
Wm = span{^^ = x; A^), 1 ^ m ^ M}. These sets are built from a subset S of V, in which the first 

sample point is picked (assuming that 0). The basis functions qm are based on and establishing 
of coefficients requires, online, the solution of a M x M system, ensuring that the EIM approximation 
wm{u, x; /x) is exact on a set of interpolation points {ti}i. 


5'i={Ri}, = w(?x(/Xi),x;/Xi), Wi = span{^i} ti = arg sup ||i(x)|, qi 


^ijx) 


( 6 ) 


Eor M p 2, the sample points fij^ are determined from a Greedy algorithm : 


Am = max^gH inRgvvM-i lk(-j •; m) - ^lli”(n) _ 

= w Pm) , Rm =-S'm-i U {/Xj^^}, Wm = Wm-i ® span{$,M} 


The system ensuring the exactness of wm-i at gives wm-i{u{Pm), x\ Pm)-, leading to the residual 

tm defined as xm^x) = w{u{Pm),x; Pm) — wm-i(u{Pm), x; Pm) on which interpolation points {tM}M ^2 and 
basis functions {( 7 m}msS 2 are computed. The next Im and Qm are given by 


Im = argsnp\rM{x)\i QmIx) 


rnjx) 

rMitn) 


( 8 ) 
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We are now ready to apply EIM on the Jacobian j the affine decomposition reads 

Qeim Q.ff 

v) + M)j‘(u,v) Vu,vGX, V/j. € V, 

q—1 m—1 1=1 

' -«-' '---^' 

non-affine part of j affine part of j 

(9) 

with (resp. the number of non-affinely (resp. affinely) parametrized terms of j, and (resp. 

j^) obtained from (resp. qi^m)- Similar decomposition is obtained for the residual r. 

2.2. Reduced basis method 

We now turn to the RB methodology. We introduce Sm = {Mi, • • • ,Mjv} with N « JV, and we define 
the set of solutions S'^ = which are orthonormalized with respect to the <, >x inner product to 

provide Wn = 

The reduced basis approximation n_/v(M) € Wn is expressed as a linear combination of Wn elements 

N 

UNifJ') = '^UN,i{f^)ii such that r{uN,v; fi; {w^{uNJx; = 0 \/v € Wn (10) 

i=l 

The offline/online strategy used to build efficiently the reduced basis approximation is based on the affine 
decomposition ([I]),®. Choosing as test functions and using (1101) we have for 1 ^ Z ^ TV 

N 

- ^UNj) = {w^{^Un , X] ( 11 ) 

1=1 

The offline step consists in the pre-computation of the terms and r^(^;),and 

is done offline once thanks to the affine decomposition. To evaluate online for any given fi, the 

coefficients Jj fi), 0* (^wat, /.t) and 7/„j(^Uiv, n), O^un^ /^) are computed in order to recover the N x N 

system giving coefficients and consequently ^+^M 7 v(/.t). 


3. A Simultaneous EIM-RB method 


The standard methodology described previously briefly requires the use of EIM prior to the RB methodology 
for each affine parameter dependent function Wi O- These EIM approximations are built up front allowing to 
then write the affine decomposition ([^l). However the greedy algorithm ([71) of the EIM offline step requires the 
computation of the solution rtjy(/x) for all points of S € T*, i.e in the step (|7|). 

The SER methodology proposes to reduce the computational cost by simply using the readily available re¬ 
duced basis approximation based on the previous EIM step and build simultaneously EIM and RB. In particular 
the expensive step ([7|) will now use solely RB approximations. 

The initialization of EIM offline stage ([Hj) doesn’t change since no reduced approximation UAr(/x) is available 
yet. The resulting rough EIM approximations Wi^i are used to compute a first affine decomposition from which 
can be built, leading to a first reduced approximation ui{fi). After this initialization step, EIM and RB 
approximation spaces are enriched alternatively. Each new EIM basis function qm ■ /*/ is built from reduced 
basis approximation UM-iip) obtained at previous iteration M — 1 (ii) complete the EIM approximation to 
then build the affine decomposition for To summarize, the EIM steps modified by SER now read 

AzM = «'’ 5 max inf ||■u;*(^^M-l(A^); •; m) “ ^llL“(n) (12) 

/xe=. zGWm-1 


^i,M = (wM-l(Mi,M);®;Mi.M) > Si^M = Wi^M = Wi^M-i®span{^i M} (13) 

The residual ri^M{x) = Wi(itM-i(/i 7 M): M^.m) “ Wi,M-iiuM-i{Pi,M)^x; then gives U^m and qi^M 


U,M = arg sup \ri^M{.x)\, 


ri,M{x) 

^i.M(ii,M) 


(14) 


The affine decomposition is then updated with the new EIM approximation Wi^M, to compute the next 
reduced basis (i.e) thanks to (fTTl) to then enrich Wn space. 
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Remark 1 (RB updates). The EIM approximations are changing during the build ofWjsr space. Then, the opera¬ 
tor - and the solved problem - evolves at each step, which can mildly deteriorate the approximation. Recomputing 
Wn elements using finite element solves based on the current EIM approximation may be considered. The offline 
precomputations need however to be updated for all elements of Wn for each update of the affine decomposition. 


Remark 2 (EIM updates). In the SER methodology only one finite element solve of ([T]) is required for the 
initialization afterwards the EIM and RB approximations are updated alternatively. We may in fact update the 
RB approximation every r EIM steps: if r = M we recover the standard method, if r = 1 we recover the SER 
method and if 1 < r < M we have an intermediairy method which requires finite element solves for the intensive 
EIM step © before the first RB update. Other alternatives are readily available and some will be discussed in 
a future publication. 


4. Preliminary results 

We now turn to numerical experiments of the SER method compared to the standard one based on the 
benchmark problem introduced in [^. It reads, find u such that 

gA‘2« _ I 

— Au + Pi- = 100sin(27rx) sin(27ry) in fl =]0,1[^ and fx = (pi, P 2 ) & T) = [0.01,10]^ (15) 

T2 

and we are interested in the output s, the average of the solution u over 12. Thanks to the function g{u, x; fx) = 
we are in the setting of the SER methodology for non-linear problems. Eollowing the standard 

M 

methodology we develop gM = lAlmix) of g thanks to the Empirical Interpolation Method (see 

2—1 

section IQ) . We shall use the absolute errors on the solution and the output, defined as follows 

=11 11^2 elpjv =1 - ■Sat I (16) 

where -jg is the finite element solution/output of the initial problem ([2]), -(v the reduced basis solution/output 
while r is the frequency at which the EIM are updated, see remark [21 

We first display the maximum of the absolute errors (ITBl) using the standard method r = M in table [Tal which 
reproduces the results in [^. Tables [T^ and (Tdl correspond to the SER methodology (r = 1) with N = M and 
investigate the influency of Wn recomputation in SER method, see remarkjTJ Table [Tbl is an intermediary stage 
such that 1 < r < M where the EIM basis functions are built by groups of r (r = 5). In this last case, since we 
do not have access yet to the RB, the first group of r EIM basis is built from finite element approximations and 
reduced basis approximations are used afterwards. This explains, for example, that (iV = 4, M = 5) is similar 
to the first one in table [Tal 

The maximum of the errors observed with SER method (r = 1) is slightly higher than with the case r = 5 
which itself displays resultats slightly higher than the standard method r = M. This behavior is expected. 
The results show the pertinence of the SER method and that we can expect good results within a reasonnable 
computational budget, since the number of finite element approximations can then be reduced to -I- 1 (table 

(Ml). 


Conclusion 

We have now an algorithm and variants that allow an efficient use of EIM and RB in the context of non¬ 
linear non-affine partial differential equations. We have already deployed SER in real applications and the 
initial results are very promising showing the same behavior as the benchmark problem of this paper. This will 
be reported as well as an analysis of SER in a subsequent paper with applications to non-linear multiphysic 
problems. 
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N 

M 

/ u.M \ ( s.M \ 

max(e^^^) max(e^_^) 

N 

M 

/ \ 
max(e^_^) 

/ s,5 \ 

maxiefjN) 

4 

5 


7.38e-3 

5.75e-3 

4 

5 

8.21e-3 

6.31e-3 

8 

10 


l.Ole-3 

2.34e-4 

8 

10 

4.48e-3 

6.18e-3 

12 

15 


1.49e-4 

3.09e-5 

12 

15 

2.69e-4 

2.36e-4 

16 

20 


2.21e-5 

1.25e-5 

16 

20 

1.48e-4 

9.31e-5 

20 

25 


5.88e-6 

2.82e-6 

20 

25 

2.60e-5 

1.46e-5 




II 




(b) r = 5 



N 

M 

/ U,1 \ 

max(e^_^) 

max(e^^^^) 

N 

M 

max(e^^^) 

max(e^^^^) 


5 

5 

9.98e-3 

7.77e-3 

5 

5 

1.30e-2 

1.02e-2 


10 

10 

2.32e-3 

1.86e-3 

10 

10 

2.20e-3 

1.50e-3 


15 

15 

4.61e-4 

3.75e-4 

15 

15 

4.83e-4 

4.05e-4 


20 

20 

2.48e-4 

2.02e-4 

20 

20 

2.42e-4 

1.98e-4 


25 

25 

3.51e-5 

2.33e-5 

25 

25 

1.50e-5 

1.24e-5 


(c) r = 1 (d) r = 1 

{Wn recomputed) (W^jv not recomputed) 


Table 1: SER algorithm: Maximum absolute errors on solution u and on output s 
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